Pygmy dipole response of proton rich argon beyond the random phase approximation 
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The occurrence of a pygmy dipole resonance in proton ricli Ar is studied using tlie unitary correlator 
operator metliod interaction V ucom, based on Argonne V18. Predictions from the random phase approximation 
(RPA) and the shell model in a no-core basis are compared. It is found that the inclusion of configuration 
mixing up to two-particle-two-holes broadens the pygmy strength slightly and reduces sensibly its strength, as 
compared to the RPA predictions. For ^^Ar a clear peak associated with a pygmy resonance is found. For ^'*Ar, 
the pygmy states are obtained close to the giant dipole resonance and mix with it. 

PACS numbers: 21.10Gv, 24.30.Gd, 24.10.Cn, 21.60.Cs, 21.60.Ev, 



I. INTRODUCTION 



The recent advances in the experimental techniques for ra- 
dioactive beams have fueled several studies of the properties 
of exotic nulcei, away from the line of beta stability. One 
of the most interesting results is the discovery of low-lying 
dipole strength in neutron rich isotopes which is interpreted as 
a pygmy dipole resonance (PDR). This new excitation mode is 
then explained as the resonant oscillation of the weakly bound 
neutron skin against the isospin saturated proton-neutron core. 
Typically, one observes in nuclei with a neutron excess, N>Z, 
a concentration of electric dipole states close to the parti- 
cle emission treshold. These carry a small fraction of the 
Thomas-Reiche-Kuhn (TRK) sum rule, which increase with 
the charge asymmetry of the nucleus. The PDR is also inter- 
esting because of its astrophysical implications. Although it 
carries only a small fraction of the total dipole strength, the oc- 
currence of increased strength at the particle separation tresh- 
old can enhance radiative capture cross sections 1 1, 2], which 
can have strong effects on the r-process nucleosynthesis and 
on the abundance distribution of elements. In addition, the 
thickness of the nuclear skin is directly related to the density 
dependence of the symmetry energy. Recently, data on the 
PDR resonance have been used to constrain models for the 
symmetry energy 

The onset of low-lying El strength has been reported even 
in stable nuclei with moderate proton-neutron asymmetry 
such as 44.48(3^ and ^o^Pb |4 H B; see Ref. |7] for a re- 
view of high precision photon scattering experiments. For 
unstable nuclei with large neutron excess, a sizable frac- 
tion of low-lying El strength was obseved in 20,22q |gj ^j^^ 
I30,i32jgjj Several theoretical models have been employed 
to study the nature of this low-energy dipole strength. Recent 
works employed the Skyrme Hartree-Fock (HP) plus quasi- 
particle random phase approximation (RPA) with phonon cou- 
plings ifTolfTlll . the quasiparticle phonon model 1 12, IslfTill . 
and the relativistic quasiparticle RPA (RQRPA) |15, In 
a recent work, the PDR was analysed including effects of 
particle-vibration coupling on top of the RPA approach UtIi . 
The inclusion of low-lying phonons mainly correspond to 
consider explicit admixtures of two-particle-two-hole (2p2h) 



states and increase the fragmentation of the dipole distribu- 
tion. In neutron rich isotopes it was found that this effect can 
generate a shift in the position of the PDR but otherwise does 
not noticeably change its characteristics. 

On the proton rich side, nuclei with an excess of protons 
over neutrons are found only for Z<50. Due to the Coulomb 
repulsion, the proton drip line is much closer to the fi sta- 
bility line and proton skins are possible only for the lightest 
isotopes. For these elements the multipole response is gen- 
erally less collective. Although these facts seem to disfavor 
the existence of a proton PDR, a recent calculation suggested 
that this mode can actually be observed in medium-mass nu- 
clei lfl6ll . Using RQRPA calculations low lying pygmy states 
were obtained when approaching the proton drip line along 
the Ar isotopes and the N=20 isotones chains. However, no 
study of correlations beyond the RPA, such as in Ref. iITtIi . 
has been reported to date for the proton PDR. In this work 
we will consider the proton rich isotopes ^^.34^ ^j^^j compare 
predictions for the dipole strength as obtained in the RPA with 
those from shell model (SM) studies in a no-core cofiguration 
space. These approaches include complementary correlation 
effects since rather different portions of the Hilbert space are 
probed. Still it is preferable to perform such a comparison 
based on the same Hamiltonian. 

In the following, the unitary correlator operator method 
(UCOM) [18] will be employed to regularize the strong core 
of the realistic Argonne V18 potential IH]. The UCOM oper- 
ator generates a unitary transformation within the many-body 
Hilbert space in which a weakly correlated wave function is 
mapped into one where strong short-range repulsion and ten- 
sor effects are explicitly manifest. If one applies this transfor- 
mation to the Hamiltonian, rather than the wave function, it 
finds an expansion of the effective interaction into many-body 
terms. The Vucom force is obtained by truncating the expan- 
sion at the two-body level. The result is an effective force 
that tames the strong short-range and tensor components of 
the original force and is therefore applicable to smaller basis 
sets. The subtle cancellation between large nuclear and ki- 
netic contributions to the total energy is also accounted for in 
this expansion. Since the coiTelator operator is chosen in such 
a way that nucleon-nucleon phase shift are preserved, Vucom 
can also be regarded as a realistic two-nucleon force on its 
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own that is applicable, however, only to medium/low-energy 
processes. Therefore, the UCOM method provides an inter- 
action independent of the employed model space and can be 
applied meaningfully to both the RPA and the SM methods. 
At the same time it includes the correlation effects proper of 
modern high precision realistic interactions. 

Details of the formalism employed in the calculations are 
discussed in Sec. The RPA and the shell model results 
are compared in Sec. Before presenting these results, we 
will discuss the dependence of our calculations on the param- 
eters of the oscillator space, and use the Hartree-Fock (HF) 
approach to construct a basis to be used in the shell model 
studies, in Sec. lIII Al Conclusions are drawn in Sec. lIVI 

11. FORMALISM 

Consistent RPA calculations were performed to guarantee 
the exact separation of the spurious center-of-mass motion. 
This means that our single particle basis was obtained by solv- 
ing the HF problem and the same interaction was used for both 
HF and RPA. To do this, one must take some care in treating 
the (i3/2 orbit which contains only two protons in the mean- 
field picture. In order to maintain consistency and to assure 
the separation of spurious states one must: (1) force uniform 
occupation for protons in this orbit when solving the HF equa- 
tions, (2) account for its partial occupation when solving the 
RPA problem, and (3) allowing for excitations of protons both 
from the ^3/2 to higher orbits and from lower orbits into the 
(half empty) level. This was done by considering the lat- 
ter orbit both as a particle and as a hole state, with the respec- 
tive depletion or occupation probabilities. Indicating the HF 
basis with a s nJajaTa ij is the isospin), the RPA eigenvalue 
equations for the amplitudes X^^^^ and are 

( RJ* AJ* )( F^'" ) " '^''( -\ )( y^'" ) ' 

where 

Kb,cd = na hh H^^i^ j^ nd n, , (2) 

/ is the angular momentum of the excited state, is the HF 
single particle energy of orbit a, and //-^^ the Pandya trans- 
form of the residual two-body interaction. The numbers («a)^ 
define the occupation of orbit a in the HF wave function, while 
(ha)^ measures the unoccupied space. Hence, for the Ar iso- 
topes we have 

II for a fully occupied 
for fl empty . (3) 

1/2 for a - K0d2/2 

In Eqs. ([TJ and (|2|, the products nafib are restricted to Iplh 
configurations only, except for the partially occupied orbit. 
It must be stressed that the above prescription arise natu- 
rally when deriving the HF-hRPA scheme from propagator 
(or Green's function) theory ll20l I21II : following Baym and 



Kadanofi^ ll22ll23ll . conserving RPA equations are derived con- 
sistently from the HF self-energy and propagator. The HF 
propagator, however, contains both particle and hole poles for 
each partially occupied orbit, such as the ii3/2 in our case. We 
checked that neglecting proton excitations from lower orbits 
to the d3/2 level no longer allows for the exact separation of the 
spurious center-of-mass mode (although the breaking is small 
and has negligible effects on the remaining El strength). The 
RPA results reported in this paper are obtained in the fully 
conserving approach. In all cases we consider the isovector 
dipole operator 

p=l 11=1 

which is corrected for the center-of-mass displacement. In 
all presentations below, the calculated B(E1) strength to each 
final state /, 

Bfi(El) = {fWm , (5) 

was folded with a Lorentzian of width r= 1 Me V, 

R(E^) = y BfiiEl) — ^ . (6) 

^ ^ (£x - w/)2 + r2/4 

For all the calculations we employ the Vucom interaction 
and subtract the kinetic energy of the center-of-mass motion. 
The intrinsic Hamiltonian is then written in the form of a two- 
body interaction, 

Hin, = r + V-T,„ = y i^^lZllL + y I . (7) 

In all cases, a UCOM operator was used which corresponds 
to a correlation volume of I§ = 0.09 fm^ for the tensor force. 
This correlator was tuned to reproduce the binding energies 
of "^He and ^H without the need of additional three nucleon 
forces ll24ll . Subsequent calculations sugggested that binding 
energies are reproduced, in perturbation theory, throughout 
the nuclear chart without additional corrections to this UCOM 
operator 1*2511261] ■ 

An harmonic oscillator basis with length bno-^-^ fm was 
used for all the RPA calculations, except for ^^Ar when in- 
dicated. This basis was truncated according to the number 
of major shells 2n + I < N„,ax (with A^„,aj— 0,1...). Calcula- 
tions based on Vucom generally converge in model spaces of 
10 to 20 shells if an harmonic oscillator basis is used. Stan- 
dard applications of the no-core shell model (NCSM) ll27ll are 
based on using harmonic oscillator bases and truncating the 
model space in terms of the maximum number of cross shell 
excitations. The reason for doing so is to allow for an ex- 
act separation of the center-of-mass motion. As it will be 
shown in Sec. IIII Al a very large number of oscillator func- 
tions would be required for converging the excitation energy 
of the pygmy state. Obviously, this is beyond the capability of 
present day computers and one needs to resort to a more real- 
istic single particle basis. Below, we will expand the Hamilto- 
nian, Eq. (|7]), over the lowest HF orbits that are obtained while 
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FIG. 1: (Color online) Isovector dipole strength of '"Ar obtained 
with HF+RPA, for different sizes of the models space (N„,a^-) and 
harmonic oscillator lengths (buo)- The dashed line represents the 
contribution from protons only [first term in Eq. Q]. 



solving the corresponding HF+RPA problem. Since we will 
only consider model spaces in which all nucleons are active, 
we refer to this scheme as NCSM. The secular problem was 
diagonalized using the shell model code ANTOINE [Ulli]. 
The dipole response was derived as usual by first applying the 
operator Q to the F=0^ ground state of ^^■^'^Ar, and then us- 
ing the resulting 1" wave function as a pivot for successive 
Lanczos iterations. 



III. RESULTS 

A. Dependence on the model space 

Figure [T] shows the dipole strength distribution, Eq. (|6]), 
obtained with RPA for different model spaces and ocillator 
lengths. In all cases one can identify a peak on the low-energy 
tail of the isovector giant dipole resonance (IVGDR) which 
receives strength almost exclusively from proton excitations, 
that is from the first term on the rh.s. of Eq. As it will be 
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TABLE I: Total El strength (in e'fm") for '"Ar obtained with 
HF+RPA for energies below 12 MeV and 40 MeV. The centroids 
of these distributions are also reported (in MeV). 



discussed below, this represents the PDR. Unfortunately, very 
large configuration spaces are needed to converge the pygmy 
peak at low energies. This behavior is an artifact of employing 
an harmonic oscillator basis, which does not describe properly 
loosely bound states. The dominant contributions to the PDR 
come from the collective excitation of protons to states just 
above the continuum threshold. Hence, rather diffuse wave 
functions may be needed to describe this strength properly. 
This implies large harmonic oscillator bases. An inspection 
of the HE single particle energies of protons show that the 
orbits in the sd shell are well converged already for N,„ax=5, 
while the pf orbits are found above the separation energy and 
keep changing when increasing the model space. 

In order to assess the sensitivity of our results to the choice 
of the model space, RPA calculations were performed up to 
Mnii.i=23 (24 major oscillator shells). For the largest model 
spaces, the HE ground state properties are independent of the 
oscillator length used. Hence, we have also preformed cal- 
culations with a larger oscillator length, bHo-2-6 fm, to fa- 
cilitates the description of states in the continuum. The cen- 
troid and total strength of the resulting El distributions are 
reported in Tab. U for energies up to 12 MeV and 40 MeV. It 
must be stressed that the pygmy behavior is reproduced by 
RPA theory even for the smallest model spaces considered 
here, and its existence is therefore a stable prediction. How- 
ever, no strength is obtained in the low-energy region when a 
small number of oscillator functions is employed. Increasing 
the space, the pygmy peak is lowered and it eventually stabi- 
lizes at around 9 MeV. Tab. |T] shows that 16 major oscillator 
shells are required to converge the centroid of the IVGDR and 
the summed strength in the region below 12 MeV. Increas- 
ing the space further, small variations of the results still oc- 
cur which are a result of discretizing the continuum with an 
increasing number of single particle orbits. The limit of an 
infinite space would be reached by solving directly the con- 
tinuum RPA equations. However, we do not expect that this 
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[MeV] 

FIG. 2: (Color online) Dependence of the HF+RPA strength distri- 
bution on the truncation of the single particle basis. The same HF 
basis is employed in both calculations, and was derived in a ^,,,01= 15 
harmonic oscillator space. Full line: all the HF states are used to 
solve the RPA equations. Dashed line; only the 6 lowest HF shells 
were retained. 



would lead to sensible deviations for the integral quantities of 
Tab.Hl 

Although very large oscillator bases are needed to generate 
the proper HF wave functions, only the orbits near the Fermi 
surface are actually relevant for constructing the pygmy res- 
onance. This is depicted in Fig. |2] were the distribution ob- 
tained in the full A^moA—lS model space is compared to the 
one obtained using the same HF basis but restricted to the 
6 lowest shells. Solving the RPA equations in this basis is 
sufficient to reproduce the original results up to excitation en- 
ergy of ~ 23 MeV . The summed strength and centroid below 
40 MeV, obtained in this way, are 22.0 MeV and 7.97 e^fm^, 
which underestimates the total strength by only a percent due 
to truncation. Thus these quantities remain largely unaffected. 
The truncation of the HF basis reproduces the pygmy results 
obtained with the larger harmonic oscillator space. It is also 
small enough to allow for no-core shell model calculations, 
which are presented in Sec. lIIICl 

We note that not only the matrix elements of the Hamilto- 
nian, Eq. Q were expanded in the new HF basis, but also the 
matrix elements of the dipole operator Q. Employing the cor- 
rect one-body matrix elements of this operator was found to be 
crucial in order to reproduce the correct transition strengths. 



B. RPA results for ^^Ar 

From Tab. U we deduce that with, the present choice 
of the correlator (Ii>=0.09fm^) and at the RPA level, the 
VucoM interaction predicts an isovector dipole strength of 
8.15±0.10 fm^e^ within the first 40 MeV of excitation energy. 
The centroid of this distribution is E/i/gdr=22.15+0.15 MeV. 
This is close to the ernpirical estimate Epeak-^^-'2/A^^^ + 
20.6/A'/*= 21.4 MeV Q. However, one should keep in 
mind that this formula was derived from data on stable nu- 
clei with masses A>50, and it is just indicative in this case. 
In the energy region up to 40 MeV the RPA calculation ex- 
haust about 153% of the TRK sum rule. This enhancement 
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FIG. 3: (Color online) Proton and neutron transition densities for 
^^Ar obtained from RPA theory (A'„,a_v=15). Three states assocated 
with the proton pygmy resonance (at 8.45, 8.60 and 10.94 MeV) are 
compared to the IVGDR eigenstates at 16.75 and 21.16 MeV . 



is expected for realistic interactions and it can be traced to 
the strong tensor component of the Vucom force ll3lll . A 
similar enhancement has also been reported applying Vucom 
to "^He (3^. Examples for the transition densities to states 
around 20 MeV are shown in Fig. [3] These describe an out 
of phase oscillation of protons and neutrons and confirm the 
IVGDR nature of these excitations. 

More interesting is the pygmy peak at the lower end of 
the GDR tail. As already noted, this receives El strength al- 
most exclusively from protons — the first term on the r.h.s. of 
Eq.© — as shown by the dashed lines in Fig.[T] The transi- 
tion densities for the RPA states at 8.45, 8.60 and 10.94 MeV 
are reported in Fig. [3] and show the typical behavior of the 
PDR: proton and neutrons move in phase in the nuclear in- 
terior while only protons are excited at the surface and ex- 
tend ouside the nuclear core. The total strength found below 
12 MeV of excitation energy is 1.0+0.15 e^fm^, for M„aA- >15. 
No strength is seen in this region if only 6 oscillator shells are 
retained. In this case the PDR is obtained at larger energies but 
it carries a similar strength (0.84 fm^e^ up to £.1 = 15 MeV). 



C. NCSM results for ^^Ar 

No-core shell model calculations were done in order to 
investigate the effect of coiTelations beyond RPA on the 
dipole distribution. As it is apparent from Fig. [T] an oscillator 
basis A^,„,n = 15, corresponding to a 13/iw model space, would 
be required to converge the energy of the PDR. The number 
of active single particle orbits can be reduced by employing 
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FIG. 4: (Color online) Top panels: isovector spectral strength ob- 
tained with the shell model for dipole excitations to the Iplh and 
2p2h configuration spaces. The first six HF orbits were included in 
the basis. Bottom panel: ITUi) result, includig excitations up to the 
pf shell only. 



the truncated HF basis discussed in Sec. IIII Al The proton 
HF orbits were used to expand the Vucom interaction and the 
dipole operator. The explicit configurations employed were 

(0il/20p3/20pi/20fl'5/2)"*""'(L?I/20fl'3/20/7/2 " ' " 0/zii/20/Z9/2)"" 

for neutrons and (0.si/20p3/20pi/20t/5/2lii/20t/3/2)'^""'' 
(O/7/2O/5/2IP3/2 •• -0/1 1 1/20/19/2)"" for protons. The model 
space was truncated in terms of the total number of particle- 
hole excitations (ny + n„ < n). This maintains all the particles 
active and allows for excitations from the core. The 0^ 
ground state was first derived using a «-particle-n-hole 
(npnh) model space. Once this is done, the dipole operator, 
Eq. (|4]i, connects this to states in a space containing up to 
(n-i-l)p-(nH-l)h. By performing Lanczos diagonalization of 
the 1" final states in the larger configuration space one is 
guarantied to exhaust the total dipole strength. Note that 
for a OpOh space, the neutrons are constrained to the lowest 
HF orbits while some mixing is possible for protons, that 
can be excited from the lower shells into the half empty ^3/2 
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TABLE II: Total dipole strength (e-fm-) and its centroind (MeV) 
in the region of the Pygmy resonace, as obtained in RPA and shell 
model with Iplh and 2p2h configurations. The dipole dirtsibution 
was integrated up to 12 MeV of excitation energy for ''-Ar and 
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orbit. This already leads to a configuration space beyond the 
corresponding RPA theory. 

The upper panels of Fig. |4] show the strength distributions 
obtained from the OpOh- Iplh and lplh-2p2h configuration 
spaces. The same structure of a IVGDR resonance and pygmy 
peak is found, like for the RPA (Fig.|2|, but both the IVGDR 
distribution and PDR are broadenend by shell model correla- 
tions. For the 2p2h model space, the IVGDR is centered at 
21.66 MeV. Integrating up to 40 MeV excitation energy, one 
finds that the total strength is lowered by configuration mix- 
ing and exhausts 137% of the TRK sum rule. The response in 
the pygmy region is compared to the RPA approximation in 
Tab. Shell model correlations do not significantly alter the 
position of the PDR. On the other hand, its strength is reduced 
by a factor of two. 

The configuration spaces used above restrict most of the 
neutrons in the (^5/2 orbit and below. One may question 
whether neutron correlations in the sd shell affect the proper- 
ties of the PDR. To check this we performed a NCSM calcu- 
lation in the Iho) model space, that is opening the sd shell and 
allowing for one excitation across major shells. The strength 
distribution, plotted in the lower panel of Fig. |4l shows no 
qualitative distortion of the pygmy peak. Hence, the ap- 
proximation made restricting neutrons in the OpOh configua- 
tion should not affects our conclusions noticeably. The PDR 
strength is also reduced here to about 60% of the RPA value. 
However, at Itio) only excitation to the pf shell are possible. 
Since this corresponds to a more severe truncation of the HF 
basis, we do not attempt to extract quantitative information on 
the PDR from this calculation. 



D. Results for ^'^Ar 

A similar analysis was carried out for ^^Ar. The strength 
predicted from RPA theory is depicted in Fig.|5]for A^mfl;t=15. 
The centroid of the dipole distribution in the interval 
[0,40] MeV is obtained as 23.3 MeV. Summing this strength 
within the same interval gives 8.69 e-fm^, which corresponds 
to 161% of the TRK sum rule. These predictions do not vary 
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FIG. 5: (Color online) Isovector spectral strength for ^*Ar obtained 
with HF+RPA (top) and the shell model up to Iplh and 2p2h (cen- 
ter and bottom). The full Ni„„^=l5 harmonic oscillator space was 
employed in the RPA calculations. The same HP basis derived for 
the RPA calculation is used for the shell model, but truncated to the 
lowest 6 shells. 



by increasing the model space and the oscillator parameter. A 
similar analysis as done for ^^Ar, suggests an uncertainty of 
+0.8 MeV and +0.2e^fm^ for these quantities. 

A pygmy resonance is also found for this isotope but closer 
to the centroid of the IVGDR. The dashed line in Fig.|5]shows 
that sizable proton contribution to the dipole response strength 
is found up to about 14 MeV. An analysis of the transition den- 
sities (Fig.|6]) confirms that these states have the character of 
a pygmy mode. Instead, the RPA eigenstates found just a few 
MeVs above contribute to the lower tail of the IVGDR. To es- 
timate the total strength associated with the pygmy resonance, 
we integrated R(Ej,) up to 14 MeV and found 0.80 e^fm^. 

No-core shell model calculations were performed in a 
model model space analogous to the ■'^Ar case, obtained by 
solving the HF equation for ^"^Ar in the same harmonic oscil- 
lator space as used for RPA. This basis was then truncated to 6 
major shells and, at the OpOh level, neutrons were constrained 
in the 1^5/2 and si/2 orbits. Also in this case we found a slight 
reduction of the total strength w.r.t. to the RPA result, yeld- 
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FIG. 6: (Color online) RPA transitions densities for the exitations to 
three pygmy states of ^^Ar (at 10.45, 12.18 and 12.71 MeV) and to a 
IVGDR state at 18.44 MeV. 



ing 8.21 e^fm^ for the lplh-2p2h configurations space. In this 
case the centroid of the IVGDR was lowered to 22.0 MeV and 
the energy weighted sum rule reduces to 143% of the TRK 
value, in the first 40 MeV. The strength resulting from the 
shell model calculation is plotted in the lower panels of Fig.|5] 
and shows a broad distribution of the El strength. Since for 
this nucleus the pygmy peak is obtained close to the giant res- 
onance, the two are strongly mixed and one finds just a long 
tail spreading in the low-energy region. From the figure it is 
not possible to clearly identify a separate peak corresponding 
to the pygmy resonance. Table |ll] lists the summed strength 
below 14 MeV, where one expects to find most of the pygmy 
states. Contrary to the ■'^Ar case, the effects of configuration 
mixing on the low-energy dipole strength is moderate. A re- 
duction of about 20% w.r.t to the RPA result was found. 



IV. CONCLUSIONS AND DISCUSSION 

The properties of the pygmy dipole resonance in the proton 
rich isotopes ^^Ar and -^^Ar were investigated by comparing 
the predictions of the RPA and the shell model in a no-core 
configuration space. The Vucom interaction was employed 
in all cases as a low-energy realistic Hamiltonian. In this ap- 
proach the nuclear force is tamed to account for the effects on 
short range coiTelations and can be directly applied in large 
scale RPA and shell model calculations. In general a large 
number of harmonic oscillator shells is still required to con- 
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verge the distribution of dipole strength in the low-energy re- 
gion. However, it was seen that the number of active orbitals 
can be significantly reduced using the HF single particle basis. 
Therefore, comparisons could be done employing the same 
Hartree-Fock basis for the RPA and shell model. 

For ^^Ar, enhanced low-lying strength was found at ener- 
gies up to 12 MeV, which could be identified as proton PDR. 
The corresponding strength distribution is peaked at about 
9 MeV, considerably above the experimental proton separa- 
tion energy of 2.4 MeV Issll . This should not come as a 
surprise, since oscillations of the proton skin need to over- 
come the Coulomb barrier. The RPA approximation predicts 
about 1.0 e^fm^ in this energy region. Additional correlations, 
as accounted for in shell model studies, reduce this strength. 
Furthermore the pygmy peak is slightly broadened due to a 
larger number of configurations. By explicitely including up 
to 2p2h states, the pygmy strength is halved to 0.49 e^fm^. 
The results obtained here predict a well defined pygmy peak 
separated from the IVGDR, in accordance with the previous 
RQRPA calculations of Paar et al. lfl6ll . Excitations associ- 
ated with a PDR were also found for ■'^Ar, which has only two 
more protons than neutrons. The present calculations place its 
strength at energies up to 14 MeV and close to the IVGDR. 
The total strength found up to this energy is also reduced by 
shell model calculations, as compared to RPA. However, the 
mixing with the nearby giant resonance states leads to a rather 
uniform response without a well defined pygmy peak. 

We note that the neutron PDR observed in larger nuclei is 
typically found below 10 MeV and its excitation energy low- 
ers with increasing mass. On the other hand, little is still 



known for medium-mass nuclei. Unstable oxygen isotopes, 
for example, show a behavior more similar to the NCSM re- 
sults of Fig.|5j i.e. a broad dipole response spreading contin- 
uously from the giant resonance region down to the PDR |8]. 
The proton PDR, if discovered, could also have a different 
qualitative behavior, as the Coulomb repulsion affects its ex- 
citation energy. The shell model predicts a centroid of the 
PDR of a;11.6 MeV for ■^'*Ar (see Tab. which may appear 
particularly large. However, this could be understood in view 
of the expected Coulomb effect. The large centroid might also 
be due to a shortcoming of the Vj/com interaction: calcula- 
tions based on this force have reported too small radii |25] and 
correspondingly they overestimate the energy of GDRs [34jl , 
at the RPA level. Although, recent second RPA calculations 
which include 2p2h configurations (as done here) seem to cure 
this issue |35]. The inclusion of three-nucleon forces may 
also play a role. Planned experiments, aimed at observmg 
low-energy dipole response of proton rich Ar isotopes f36ll . 
will help testing the accuracy of the present approach. If the 
pygmy states in Ar are actually observed within a few MeVs 
of the IVGDR tail, a substantial mixing should be expected. 
This would result in a broad PDR, seen as an extended low- 
energy tail of the isovector dipole distribution. 
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